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ABSTRACT 

We investigate the attenuation law seen through an interstellar medium (ISM) with 
clumpy spatial distributions of stars and dust. The dumpiness of the dust distribution 
is introduced by a multi-phase ISM model. We solve a set of radiative transfer equa- 
tions with multiple anisotropic scatterings through the clumpy ISM in a 1-D plane- 
parallel geometry by using the mega-grain approximation, in which dusty clumps are 
regarded as very large particles (i.e. mega-grains). The dumpiness of the stellar dis- 
tribution is introduced by the youngest stars embedded in the clumps. We assume 
a smooth spatial distribution for older stars. The youngest stars are surrounded by 
denser dusty gas and suffer stronger attenuation than diffuse older stars (i.e. age- 
selective attenuation) . The apparent attenuation law is a composite of the attenuation 
laws for the clumpy younger stars and for the diffuse older stars with a luminosity 
weight. In general, the stellar population dominating the luminosity changes from 
older stars to younger stars as the wavelength decreases. This makes the attenua- 
tion law steep; the composite attenuation rapidly increases from small attenuation 
for older stars at a long wavelength to large attenuation for younger stars at a short 
wavelength. The resultant attenuation law of normal disc galaxies is expected to be 
much steeper than that of starburst galaxies observed by Calzetti et al. Finally, the 
Calzetti's attenuation law is regarded as a special case with a large density in our 
framework. 

Key words: dust, extinction — galaxies: ISM — galaxies: spiral — ISM: structure 
— radiative transfer 



1 INTRODUCTION 

Since the first dust grain condensation in the universe, dust 
grains have accumulated in the interstellar medium (ISM) of 
galaxies. This internal dust attenuates the stellar light be- 
fore it escapes from galaxies. Thus, we need a correction of 
the internal dust attenuation in order to obtain the intrinsic 
flux and the spectrum of galaxies. This procedure should be 
done with an appropriate attenuation 1 law including effects 
of scatterings and the configuration of dust and stars. Gen- 
erally, it is different from the extinction laws observed in the 
Milky Way and Magellanic Clouds; an extinction law corre- 
sponds to the attenuation law in a distant uniform screen 
geometry which is not realistic at all for most galaxies. How- 
ever, there are many works in which an extinction law is used 
to correct for the internal dust attenuation. 
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1 In this paper, we use the term, "attenuation" . However, other 

terms (e.g., "e ffective extinct ion" , "obscuration") are also used in 

the literature <Calzettill200ll) . 



Observationally, only an average attenuation law of 
ultra-violet ( UV) bright starbur st galaxies has been ob- 
tained so far dCalzetti et alJll994 hereafter it is called the 
Calze tti law, see also lLertherer et alJl2002l and iBuat et alJ 
2002). This is because the UV spectra o f galaxies were rare . 
GALEX (GALaxy Evolution eXplorer; iMartin et al"ll2005l) 
will change the situation. We will have UV spectra of a lot 
of galaxies. From the data, we will obtain attenuation laws 
of various kinds of galaxies. 

Theoretically, to study the attenuation law of galaxies 
is to solve the radiative transfer in a galactic scale. There 
are several groups studying the galactic radiative transfer 
(see ICalzettil 1200 ll for a review and references therein). 
We will present a brief summary for the works about the 
attenuation l aw below. Because o f the importance of scat- 
terings (e.g., iBruzual et al.lll988l) . we restrict ourselves to 
only works which include scatterings. 

The attenuation law through a smooth dust and stel- 
lar distribu tion have been exam i ned since the early phase of 
this topic iBruzual et alJll988t iDi Bartolomeo et al.l fl99^: 
iFerrara et al.ll999l : lBaes fc Deionghel2001bl) . In general, the 
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attenuation amount is smaller than that by a distant uni- 
form screen because of the geometric and scattering effects. 
For example, let us compare a plane-parallel slab where dust 
and stars are well mixed with the foreground uniform screen 
which has the same amount of dust as the slab. Obviously, 
the optical depth from an outside observer to stars close to 
the surface of the slab is smaller than that of the foreground 
screen. This is a geometric effect. The scattering effect is 
that some scattered photons coming into an observer's line 
of sight compensate a part of the radiation absorbed and 
scattered out from the line of sight. Therefore, the geomet- 
ric and scattering effects make the medium less opaque. 

The dumpiness of the medium (i.e. dust distribution) 
is als o important as clearly shown by iNatta fc Panaeial 
Jl984T) . This also makes the medium more transpar- 
ent than a uniform screen geometry because there 
are photons e scaping from the medium through a less 
opaque path iGordon et alJ Il997l: Varosi fc Dwekl Il999t 
[ Gordon et all I200C : IWitt fc Gordon! l200fJ: iBianchi. et alJ 
120001: iPierini et alJl2004H ~ Since such effects are more effi- 
cient for shorter wavelengths, the resultant attenuation law 
usually becomes grayer th an the extinction law (distant uni- 
form screen). For example. lGordon et alJ il997t) showed that 
the attenuation law through the clumpy medium with the 
Small Magellanic Cloud (SMC) type dust (i.e. no prominent 
feature at 2175 A) can be very similar to the Calzetti law 
which is much grayer than the assumed SMC extinction law. 

We can also consider t he dumpine s s of t he stellar dis- 
tribution. As introduced bv ISilva et al.l (Il998t) . young stars 
may localize in their birth clouds (i.e. molecular clouds). Be- 
cause of a larger density of the surrounding gas, the radiation 
from young stars suffers stronger attenuation than that from 
stars distributed diffusely outside the birth clouds. With this 
hypoth esis (hereafter it is called the age-selective attenua- 
tion), |Gj^nato_et_alJ B200Ci) well reproduced the Calzetti law 
with the Milky Way (MW) type dust (i.e. with the 2175 A 
bump). For starburst galaxies, the radiation field is domi- 
nated by the young stellar radiation which suffers a feature- 
less and grayer attenua t ion la w expected through a dense 
medium. iGranato et al.l a2000f) also obtained a different at- 
tenuation law for normal galaxies w here the diffuse compo- 
nent plays a role. iTuffs et alJ <l2004l) showed that the global 
attenuation law is a composite of attenuation laws of various 
stellar components with weights of their fluxes. 

However, there has been no study of the attenuation law 
taking into account the d umpiness for both dis tributions of 
dust and stars, except for lBianchi. et ail j200(J) . Since their 
analysis is monochromatic, it is not the attenuation law in 
a strict sense. In real galaxies, the ISM is very clumpy, and 
the localization of young stars is often observed. Therefore, 
we should treat the dumpiness for both of stars and dust to 
develop the investigation of the attenuation law. The current 
paper executes it for the first time. 

This paper also discusses the relation between the 
dumpiness and the physics of the medium. In the radia- 
tive transfer through the clumpy two-phase medium, there 
are two essential parameters: the density contrast between 
clumps and the inter-clu mp medium, and the v olume filling 
fraction of clumps (e.g., IWitt fc Gordon! Il996f) . We trans- 
late these two parameters into more fundamental quanti- 
ties, the ISM thermal pressure and t he mean density, b ased 
on the multi-phase ISM model (e.g.. IpTeld et al.lll969l) : we 
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Figure 1. Phase diagram: thermal pressure— hydrogen number 
density. The dotted curves indicate the thermal equilibrium 
points in the interstellar medium of the Milky Way (top: Galac- 
tocentric radius of 3 kpc, middle: 8.5 kpc, and bottom: 15 kpc) 
reproduced from fig. 7 of Wolfire et al. (2003). The thick solid lines 
are the approximate relations of two thermally stable phases; the 
warm neutral medium (WNM, left) and the cold neutral medium 
(CNM, right). Their analytical forms are given by equations (1) 
and (2). 



try to relate the attenuation law with th e ISM physics. 
Recen tl y, such an attemp t was started bv iFischera et all 
( 2003): iFischera fc Dooital l)2005|h Since their analysis was 
only a distant foreground screen (not uniform but turbulent) 
case, the geometry considered in this paper is more suitable 
for extended sources like normal disc galaxies. 

This paper is organized as follows: in section 2, we intro- 
duce a physical model of the dumpiness as an application 
of the two-phase ISM model. Then, the radiative tr ansfer 
and the mega-grain approximation llVarosi fc Dwekl ll999T> 
are described in section 3. The main result of this paper, a 
steep attenuation law, is presented in section 4. The robust- 
ness of the result, the Calzetti law, and a simple prescription 
of the dust attenuation proposed previously are discussed in 
section 5. The final section is devoted to the conclusion of 
this paper. 



2 CLUMPINESS OF THE ISM 

The ISM in galaxies is very clumpy, and the dumpiness 
is related with the physics of the ISM . As an applica- 



tion of the multi-phase ISM model (e.g., iField et alJ 119691 : 
iMckee fc Ostrikerll977l) . a physical model of the dumpiness 
is introduced in this section. 

Assuming thermal energy and chemical equilibria in the 
ISM with a temperature less than 10 4 K, we find two ther- 
mally ^table^hAsesinthe^^ 

fe.g.. lWolfire et alJll995t iKovama fc Inutsukall2000D : one is 

the warm neutral medium (WNM) and the other is the 
cold neutral medium (CNM). Figure 1 is an example of the 
phase diagram. The dotted curves indicate the equilibrium 
po ints in environments of the Milky Way's ISM calculated 
bv lWolfire et alJ (120031) . These curves depend on the heat- 
ing and cooling functions considered; they depend on the 
metallicity, the dust amount and properties, the intensity 
of the radiation field, the flux of cosmic rays, and so on. 
However, differences are apparent mainly in the local maxi- 
mum and minimum pressures; locations of thermally stable 
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phases (especially for the WNM) are rather robust. Thus, 
we adopt the following approximate relations between the 
thermal pressure and the hydrogen number densities in the 
two phases: 



Table 1. Physical quantities determining the attenuation law. 
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nn,- 



10 4 Kcm" 3 
and 

p/k B 



1 cm" 



(WNM), 



(it&r (CNM). 



(1) 



(2) 



10 45 Kcm- 3 V10 3 cm- 
They are shown in figure 1 as thick solid lines. 

We also assume a thermal pressure equilibrium among 
phases. Although the total ISM pressure is dominated by 
the non-thermal (turbulen t motions, magnetic f ields, and 
cosmic rays) pressure (e.g.. lBoulares fc CoxlH^ . thermal 
pressures of warm and cold neutral gases and hot gas are 
simil ar observationa ll y, p/k g = 10 3-4 K cm -3 (e.g.. iMversI 
ll97SMWoffire et alJ (2003) proposed a scenario where the 
thermal energy and pressure equilibria, and then, the two- 
phase neutral medium should be established in the Milky 
Way's ISM even in a turbulent condition. Such equilibria are 
probably valid for the ISM in other normal disc galaxies, at 
least as a long term averaged property. 

In the radiative transfer through a two-phase medium 
discussed here, we regard the WNM and the CNM as the 
inter-clump medium and clumps, respectively. Here, the to- 
tal filling fraction of the WNM and the CNM is assumed 
to be un ity for simplicity althou gh it is about 0.5 observa- 
tionally (iHeiles fc Trolandll2003h . In fact, the rest ~50% of 
the volume is filled with the hot (~ 10 6 K) gas produced by 
supernovae. This hot gas determines the thermal p ressure of 
the WNM and the CNM iMckee fc Ostrikerlll977ft . 

Once an equilibrium thermal pressure in the medium 
is given, we have the densities of the CNM and the WNM 
based on equations (1) and (2). In other words, we have 
the density contrast between clumps and the inter-clump 
medium. Then, we assume a mean density of the hydrogen 
nucleus which gives the volume filling fraction of clumps; 

/cl = (nn — HH lW nm)/(nH,cnm — ^H.wnm)- 

Here, we deal with normal disc galaxies. An equilibrium 
therm al pressure of p/k-g. = 10 3 ' 5 K cm -3 dWolfire et alJ 
120031 . for the solar vicin ity) and a me an hydrogen number 
density of riH = 1 cm" " 3 JSpitzerlll97Sl) are assumed. These 
values are listed in Table 1 where other quantities assumed 
later are also summarized. The corresponding density con- 
trast and clump filling fraction are 118 and 1.85 x 10 -2 , 
respect ively. This density cont rast is similar to that consid- 
ered in lWitt fc Gordon! <|2000T) . whereas our filling fraction 
is about an order of magnitude smaller than theirs. That is, 
our mean density of t he medium is much less than that in 
IWitt fc Gordon! (l2000l) . 

Finally, clumps are assumed to be self-gravitating; the 
clump size, which is required to calculate the optical depth 
of each clump later, is given by the equilibrium pressure 
and the CNM density through the Jeans length. This is an 
assumption motivated by the following discussion where we 
regard clumps as birth clouds of young stars. Indeed, the 
molecular clouds should be in the CNM. For the equilibrium 
pressure and mean density assu med abov e, the Jeans length 
for a sphere of the CNM is ^/l5p/4nG/p c i = 10.4 pc (= 
r c i), where p is the thermal pressure, G is the gravitational 
constant, and the density inside clumps p c i = /im p nH,cnm 



Name 


Notation 


Standard value 


equilibrium thermal pressure 


pAb 


10 3 5 K cm" 3 


mean hydrogen density 


"H 


1 cm -3 


dust-to-gas mass ratio 


V 


icr 2 


half height of the gas+dust disc 




150 pc 


layering parameter of old stars 




0.5 


young stellar age criterion 


ty 


10 7 yr 


homogeneous visual optical depth 


(face-on) 


0.54 



with the mean atomic weight p = 1.4 and the proton mass 
m p . This size is similar to (but sli ghtly larger th an) a typical 

3nl ll98ll) . 



size of the molecular clouds (e.g.. lLarso: 



3 RADIATIVE TRANSFER 

3.1 Dust properties 

In this paper, we use an empirica l set o f dust optical prop- 
erties taken from lWitt fc Gordon! |2000). This data set con- 
sists of the dust opacity (kd,\), the scattering albedo (ujd,x), 
and the asymmetry parameter (gd,\) for the MW type and 
the SMC type dust. The number of wavelengths is 25 from 
0.1 [im to 3 [im. Since there is only the wavelength de- 
pendence for kd,x in the data set, we adopt the absolute 
value of the opacity at visual band as kd,v = 2.5 x 10 4 
cm 2 g" 1 (ext inction cross section per unit dust mass; e.g., 
lDraindl2003T) for both types of dust although the actual 
values for the MW dust and the SMC dust are slightly 
different. Effects of the size distribution and the compo- 
sition are already included in the data. For calculations 
of the scatterin g angle, we assume th e Henyey-Greenstein 
phase fu nction j Henve v fc Greensteinl Il94 ll) for simplicity 
although iDraind a2003l) proposed a better function. Finally, 
we assume the dust-t o-gas mass rat io in the ISM of the Milky 
Way V = 1Q- 2 (e.g.. lSpitzerlll978t) . as a standard value. 

3.2 Mega-grain approximation 

To treat dumpiness in the radiative transfer problem, we 
have to consider a spatial 3-D geometry. However, we can 
do in a spatial 1-D ge ometry w i th the mega-grain approx- 
imatio n introduced by iNeufeldl l[l^9lJjj jH^bsonfc_l^admanl 
lll993l) and further developed bv ]^^o^^^DwekT^^99l) . In 
this approximation, we regard a dusty clump as a huge parti- 
cle called a mega-grain which produces absorption and scat- 
tering effects like a normal dust grain. Practically, we replace 
usual optical properties of normal dust grains (extinction co- 
efficient, albedo, and asymmetry paramete r) with effective 
ones. We follow the formulation given by IVarosi fc Dwekl 
(1999), in which an extensive comparison between the ap- 
proximation and 3-D Monte- Carlo radiative transfer calcu- 
lations has been made in some spherical geometries and the 
validity of the approximation has been shown clearly. We 
give a summary of equations used in this paper below. We 
will omit to write the wavelength dependence explicitly. 

We assume that all clumps are an identical sphere with 
the radius r c \ and the gas density p c \ given in section 2 and 
these clumps are floating in the inter-clump medium with 
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the gas density pi cm = A tm p n H,wnm- The optical depth (ra- 
dius) of a clump relative to the inter-clump medium is 



Tci = (pd — picm)kdVr c \ ■ 



(3) 



This optical depth vanishes in the case of no density con- 
trast. Since the interaction (absorption and scattering) prob- 
ability against the parallel light by a sphere with an optical 
depth (radius) r is exactly 



P int (r) = l-^ + 



(7 



1 

2^ 



(4) 



the extinction coefficient per unit length of the medium by 
clumps (i.e. mega-grains) is 



= n c l7T7-;?i-Pmt(Tcl) = ^7^- Pint (t, 

4r c i 



(5) 



where the clump number density n c \ was replaced with 
3/ci/47rr^. We did not consider any overlap of clumps. Then, 
the effective extinction coefficient per unit length in the two- 
phase medium is defined as 



Kcff 



+ Kic 



(6) 



where Ki cm = fcdDpicm is the extinction coefficient per unit 
length of the inter-clump medium. 

The effective albedo u) e g and the effective asymmetry 
parameter g c s are also defined as 



Weff ; 

and 

Sefi = 



tteff 



Pelting "t" <?dticm 
Kcff 



(7) 



(8) 



respectively, where u c i and g c \ are the albedo and asymmetry 
parameter of a clump, respectively. IVarosi fc Dwekl l|l999Tl 
proposed an expression of uj c \ as 



U>c\ — Wd-PeBc(T c l, Wd) , 



where 



P cbc (t,uj) = 



(3/4r)P int (r) 



l-w[l- (3/4r)P int (-r)] 



(9) 



(10) 



This is the photon escape probability from a sphere in which 
isotropic sources and dust whose scattering is isotropic dis- 
tribute uniformly. In the derivation of equation (10), the 
uniform distribution of the sc attered photons is als o as- 
sumed (see Appendix C of IVarosi fc Dwekl I1999T) . Al- 
though equation (9) is an approximation, the agreement 
with results obtained b y the Monte Carlo simulation is good 
JVarosi fc DweJdll99Sh . For the clump asymmetry parame- 
ter, we use equation (43) of lVarosi fc Dwell lll999f) which is 
an analytical fitting formula of their Monte Carlo simulation 
results. 

In the case with a large fi lling fraction of clump s / c j , 
the overlap of clumps happens. IVarosi fc Dwekl lll999l) pro- 
posed the replacement of r c \ with r c i(l — , f c i) 7 t o take into 
account this overlap effect. IVarosi fc Dwekl <ll999T) also noted 
that 7 = 1 is suitable for the case of uniformly distributed 
sources which is our case. Although it is included in the com- 
putation, this effect is almost negligible for the standard case 
discussed later because / t i < 1. 



3.3 Transfer equations 

Suppose a 1-D plane-parallel geometry along z-axis. The 
mid-plane is set to be z = where we put a mirror bound- 
ary. Above the mid-plane, we put the gas+dust disc with 
a constant mean density up to a height of /id- Here, we 
assume ha = 150 pc which is the observed height of the 
cold neutral and molecul ar gas in the Milky Way (e.g., 
iBinnev fc Merrifieldl Il99sh . Although its mean density is 
constant, it has a clumpy structure described in section 2, 
and its extinction coefficient is the effective one K e g given 
by the mega-grain approximation. For simplicity, we assume 
that there is no dust outside the gas+dust disc. Then, the 
optical depth coordinate is defined as dr — — n c sdz with 
t = at z = hd' Since K e g is constant throughout the 
gas+dust disc, we have r = n^ghd at z = 0. Note that the 
effective albedo ui a s and the effective asymmetry parameter 
g e ff are also constant throughout the disc. 

A static radiative transfer equation in this geometry is 



H^^- = S(T,„)-I(r, f ,), 



(11) 



where I is the specific intensity along a ray and /1 is the 
cosine of the angle between the ray and the z-axis. The 
source function S is given by 

S(t, 11) — rj»(r)/K e B + w e ff / I(t, n')$(g e s, p, fi')dfi , (12) 



where T], is the stellar emissivity (isotropic) and $ is the 
scattering phase function. Note that u; e ff and g c s appear in 
the second term of the right hand side (i.e. scattering term). 
The boundary conditions are 



at z — (mirror boundary) and 

/(r = 0,M<0) = -^, 
M 

at z = hd, where 
I* = r/ t (z)dz . 

Ah 



(13) 
(14) 

(15) 



3.4 Emissivity distribution — Stellar dumpiness 

We consider two kinds of the emissivity source: one is em- 
bedded in clumps and the other distributes diffusely and 
smoothly. Then, we consider that the embedded source is the 
youngest stars (the age criterion is defined later) and other 
older stars distribute diffusely. This is motivated by the fact 
that stars are formed in molecular clouds (i.e. clumps) and 
diffuse as ageing. This makes the dumpiness of the emissiv- 
ity. Here we normalize the intrinsic (i.e. in the no dust case) 
emissivity of each component as 



rj t (z)dz 



3.4-1 Clumpy stellar emissivity 



(16) 



When clumps distribute uniformly in the gas+dust disc (we 
did not consider any spatial correlation of clumps in section 
2), a mean density of stars embedded in clumps is constant 
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throughout the gas+dust disc. On the other hand, there 
is no embedded star outside the disc because there is no 
clump. Thus, the intrinsic clumpy stellar emissivity is 1/2/id 
for \z\ ^ /id and zero for \z\ > hd when it is normalized 
as equation (16). A part of the radiation energy from the 
embedded stars is locally absorbed by dust in a clump. If 
stars distribute in the clump uniformly and the scattering 
by dust is isotropic, we can adopt the escape probability of 
equation (10). Therefore, the clumpy stellar emissivity input 
into the transfer equation is 



/ \ I ^esc 

= \ o 



(Tci,Wd)/2ftd (for \z\ ^ h d ) 
(for \z\ > h d ) 



(17) 



In this case, we have I, = 0. 



3.4-2 Diffuse stellar emissivity 

The diffuse stellar disc is assumed to have an exponential 
structure along z-axis with a scale height of hd/£, where £ 
is the layering parameter, i.e. the ratio of the heights of the 
dusty disc to the stellar disc. Here we set £ = 0.5 which 
means that the stellar scale height is t wice times larger than 
the h eight of the gas+dust disc (e.g., iBinnev fe Merrifieldl 
1998). From the normalization of equation (16), we have 



2h d 



(18) 



72- 



3.5 Computational scheme 

The radiative transfer equation described above was solved 
by the A-iteration s cheme with a difference formula of the 
formal solution (e.g.. lMihalas fe Weibel-Mihaladll999l) . The 
iteration was continued until the maximum fractional dif- 
ference between previous and present values of S becomes 
less than 10 -10 . To accelerate the convergence, we adopted 
the Ng's algorithm jNdll974l : Iplson et alJll986T) . Thanks to 
a modest optical depth of the problem discussed here, the 
number of iterations was always less than 10-20 times. 

As a code test, we compared the inclination depen- 
dences of the attenuation amount for several optical bands 
in a smooth plain-parallel geometry calculated by the ra- 
diative transfer code of this paper with that calculated 
by a spherical harmonic s method presented in table 2 of 
iBaes fe Deionghd (l2001aH . For the comparison , we ad justed 
the code to the set up of IBaes fe Deionghd J2001a). The 
agreement was excellent: differences were within a few % in 
magnitude. 



3.6 Transmission rate 

We first define the transmission rate for each emissivity 
source. The intensity observed by a distant observer (in a 
direction n > 0) is 



/ obs Gu>0) = /(T = 0,/x>0) + 



I' 



(19) 



On the other hand, the intrinsic (no dust case) intensity to 
the observer is 



galactic age lOGyr 

i y =10Myr . 

i y =100Myr 

£ y =lMyr 




0.1 



1.0 



X [/zm] 



Figure 2. Luminosity fraction of young stellar population as a 
function of the wavelength at the galactic age of 10 Gyr. The solid, 
dashed, and dotted curves are the cases with the age criterion 
for young stars of 10 Myr, 1 Myr, and 100 Myr, respectively. A 
constant star formation history is assumed for thick curves, and 
an exponentially decaying star formation history with 5 Gyr time- 
scale is assumed for thin curves. The difference between the two 
dashed curves are too small to see it. 



Iint(fJ- > 0) = - , 



(20) 



because of the normalization of equation (16). Thus, the 
transmission rate is 

7V > 0) = ^ g § = Hl{r = 0, n > 0) + J. . (21) 

When we consider two components of the emissivity source 
(clumpy and diffuse), the total transmission rate becomes a 
composite of transmission rates for them as 



T = fTi + (1 - f)T , 



(22) 



where / is the fraction of the emissivity embedded in clumps, 
and 

Ti: The transmission rate for the clumpy source. 
To: The transmission rate for the diffuse source. 

Hereafter, these abbreviations are used. 

Since we consider that the embedded source is the 
youngest stars and the diffuse source is older stars, Xi and T 
are transmission rates for the radiation from younger stars 
inside clumps and for the radiation from older stars outside 
clumps, respectively. The emissivity fraction / is now the 
luminosity fraction of the younger stars. We denote it as f y . 
It depends on the star formation history in addition to the 
age criterion of young stars. If the star formation rate at a 
galactic age t is denoted by *&(£), the luminosity fraction f y 
as a function of the age t is 



fy(t) = 



f™ nlt *' tl y(t-t')L(t')dt' 
J^(t-t')L(t')dt' 



(23) 



where t y is the criterion of the young stellar population, and 
L{t) is the luminosity density of the simple stellar population 
at the age t. 

Figure 2 shows the wavelength dependence of / y for 
several values of t y in the galactic age of 10 Gyr case. Al- 
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though we assumed two star formation histories, a constant 
star formation rate and an exponentially decaying star for- 
mation with 5 Gyr time-scale, the differences are very small. 
Hereafter, we discuss only the constant star formation case. 
On the other hand, the choice of t y is more important. 

The age criterion, t y is physically related with the life- 
time of the birth clouds (i.e. clumps) and/or the time- 
scale of stars moving out from the birth clouds. The life- 
time of the molecular c louds is estimated to be about 10 
Myr (Blit z~ Shulll980l) . The moving off time-scale is also 
about 10 Myr for the size of clumps considered here since 
the obser ved velocity disp ersion of y oung stars is ~ 1 km 
s" 1 (e.g-. lde Bruiinelll99gl) . Moreover. ICalzetti et all (Il994l) 
found larger attenuation against the hydrogen recombina- 
tion lines than that against the UV stellar continuum. This 
implies that ionizing stars (i.e. O-type stars) are more deeply 
embedded in dusty clouds than later-type stars. The main- 
sequence life-time of the latest O-type stars is also about 10 
Myr. Therefore, we choose 10 Myr as a typical value for t y 
JSilva et alJll998tlCharlot fc FalfeoOOl) . 



4 RESULT: A STEEP ATTENUATION LAW 

We first show a typical attenuation law of normal disc galax- 
ies with the standard values of physical quantities summa- 
rized in table 1. For comparisons with the literature, we note 
the equivalent face-on visual optical depth in the homoge- 
neous case: 




for the visual extinction cross section per unit dust mass 
k i)V = 2.5 x 10 4 cm 2 g" 1 . 

Figure 3 shows the transmission rates for the face-on 
(n = 1) case: the panel (a) is the case assumed the MW 
type dust (with the 2175 A bump) and the panel (b) is the 
case assumed the SMC type dust (without the bump). The 
solid curves are the cases of the clumpy dust distribution: 
in each panel, the upper thin curve is the transmission rate 
To and the lower thin curve is the transmission rate Ti. 
The thick solid curve in each panel is the typical case for 
normal disc galaxies, i.e. the composite transmission rate by 
equation (22) with f y displayed as the solid curve in figure 2 
(t y = 10 Myr). For a comparison, the uniform screen case 
and the uniform gas+dust disc case (i.e. no dumpiness for 
stars and dust) are shown by the dotted and dashed curves, 
respectively. 

Generally the ISM dumpiness makes the medium less 
opaque. Indeed, the transmission rates increase in order of 
the screen (dotted), the uniform disc (dashed), and the 
clumpy disc (upper thin solid, To). However, we can have 
heavier attenuation if s tars are embedded in clumps (lower 
thin solid, Ti) (see also lBianchi. et aljEciooh . 

The composite transmission rate, T\ follows To,a at long 
wavelengths where old stars outside clumps dominate the lu- 
minosity (i.e., / Yi a — 0), whereas it approaches T lt \ at short 
wavelengths where young stars inside clumps dominate the 
luminosity (i.e., f y ,\ ~ 1). This makes the attenuation law 
steep. As shown in figure 4, the expected typical attenua- 
tion law of normal disc galaxies (thick solid curve) is much 



steeper than the Calzetti law 2 (dot-dashed curve), the ex- 
tinction laws (i.e., uniform screen case, dotted curve) and 
the attenuation law for the uniform stellar and dust distri- 
bution (dashed curve). This is the main conclusion of this 
paper. 

Before moving to the examination of the parameter de- 
pendences, we note that the transmission rate slightly larger 
than unity for a very small opacity case (e.g., longer wave- 
lengths). For the case displayed in figure 3, the transmission 
rate is 1.004 at 2 [im, that is, 0.4% brighter than the in- 
trinsic luminosity. This i s not artefact but real d ue to the 
effect of scatterings (e.g.. lBaes fc Deiongheil2001bl see also 
section 5.1.7). 



5 DISCUSSION 

5.1 Robustness of the result 

A steep attenuation law obtained in the previous section is 
mainly caused by the age-selective attenuation; a composite 
of two different attenuations, namely, smaller attenuation for 
diffuse stars and larger attenuation for young stars embed- 
ded in clumps. Here, we examine which parameters affect 
on the conclusion by changing values of them from the stan- 
dard case (table 1) one by one. In figures of this section, 
we will see only one case of the MW type dust or the SMC 
type dust because the two dust cases give qualitatively very 
similar results. 

5.1.1 Young stellar age criterion 

First we change the criterion of young stars, t y . As shown 
in figure 2, it affects on the luminosity fraction of young 
stars, / y . Figure 5 shows the effect of t y on transmission 
rates and attenuation laws. For the thick solid, dotted, and 
dashed curves, t y = 10 Myr, 100 Myr, and 1 Myr are as- 
sumed, respectively. Two thin solid curves in the panel (a) 
are transmission rates To (upper) and Ti (lower). As t y in- 
creases, / y increases and approaches unity, so that the com- 
posite transmission rate decreases and approaches T\. On 
the other hand, the attenuation laws normalized at visual 
band for t y > 10 Myr are very similar. The conclusion of a 
steep attenuation law is not much affected by t y if t y > 10 
Myr. Even for the case of t y = 1 Myr, the attenuation law 
is steeper than the Calzetti's law. 

5.1.2 Layering parameter for diffuse stellar disc 

Even if we change the layering parameter for diffuse stellar 
disc £, the resultant transmission rates and attenuation laws 
are almost not changed as shown in figure 6 where three 
cases of £ = 0.5 (solid), 1.0 (dotted), and 2.0 (dashed) are 
displayed. This is because the typical disc considered here is 
not so opaque (the visual optical depth is about 0.5 in the 
uniform screen geometry). If we considered an opaque disc, 
the transmission rate (equation 21) would be determined by 
the second term which depends on the layering parameter £ 
more strongly, so that some differences would appear. 

2 We extrapolat ed the attenuation la w presented by ICalzettil 
bOOllK based on lLeitherer et alj <2002l) for A < 0.18um. 
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Figure 3. Face-on transmission rates for (a) the Milky Way type dust and (b) the Small Magellanic Cloud type dust. In each panel, 
the upper thin solid curve is the transmission rate for the source outside clumps (To), the lower thin solid curve is that for the source 
inside clumps (Ti), and the thick solid curve is a composite of the two cases with the luminosity weight f y displayed as the solid curve 
in figure 2 (the age criterion of young stars is 10 Myr). The dotted and dashed curves are the uniform screen and the uniform disc (no 
dumpiness for stars and dust) cases, respectively. The height of the dusty disc is 150 pc and the layering parameter of the diffuse stellar 
disc is 0.5. Other assumed quantities are displayed in the panels. 
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Figure 4. Face-on attenuation laws (normalized at visual band) for (a) the Milky Way type dust and (b) the Small Magellanic Cloud 
type dust. The solid curve is the typical case for normal disc galaxies denoted as the thick solid curve in figure 3. The dotted and dashed 
curves are th e uniform screen and the uni form disc (no dumpiness for stars and dust) cases, respectively. The dot-dashed curve is the 
Calzetti law iCalzettill200 li iLeitherer et al . 2002). The height of the dusty disc is 150 pc and the layering parameter of the diffuse stellar 
disc is 0.5. The criterion of the young stellar age is 10 Myr. Other assumed quantities are displayed in the panels. 



5.1.3 Inclination angle 

Figure 7 shows the composite transmission rates and at- 
tenuation laws for several cases of the inclination angle. As 
the inclination angle increases (ft decreases), the disc opti- 
cal depth increases because the disc apparent thickness has 
a dependence of 1/jJ,. The effect of the disc edge appears 
from an inclination angle larger than 89° if the ratio of the 
disc radius to the vertical thickness is 100. However, the ra- 
dial density structure may affect on a highly inclined case. 
When we restrict ourselves to a inclination angle of 45° , for 



example, the difference from the face-on case is still small, 
i.e. a steep attenuation law. 



5.1.4 Mean gas density 

Figure 8 shows the effect of the mean gas density on the 
composite transmission rates and attenuation laws. To keep 
the two-phase medium described in section 2, the hydro- 
gen number density is limited in the range of rm.wnm ^ 
n-R ^ HH,cnm- As the hydrogen number density increases, 
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Figure 5. Dependence on the criterion of the young stellar age 
(ty). (a) transmission rates and (b) attenuation laws normal- 
ized at visual band. In both panels, thick solid, dotted, and 
dashed curves are the composite transmission rates/attenuation 
laws with t y = 10 Myr, 100 Myr, and 1 Myr, respectively. Thin 
solid curves in panel (a) are transmission rates To (upper) and Ti 
(lower). The Milky Way type dust is assumed. 



the transmission rates decreases and the attenuation laws 
becomes grayer because the disc optical depth increases. 3 
For normal disc galaxies, a mea n hydrogen nu mber density 
is typically around 1 cm -3 (e.g.. ISpitzerlll97ci) in which we 
expect a steep attenuation law. On the other hand, a larger 
mean density is expected for actively star-for ming galaxies, 
e.g., U V bright starburst galaxies observed bv lCalzetti et all 
( 1994). This means that the Calzetti law can be reproduced 
by a high density condition. Indeed, this is found in section 
5.2. 



5.1.5 Equilibrium thermal pressure 

When we change the equilibrium thermal pressure, we 
should direct our attention to the local minimum and max- 
imum pressures. Since we stand on a two-phase ISM model 
described in section 2, we restrict ourselves within the range 
between p/ks = 10 3 K cm -3 and 10 4 K cm -3 expected that 
a two-phase medium is established in the environment o f 
the Milky Way's ISM (figure 1, see also lWolfire et al.ll2003D . 
We show the dependence on the thermal pressure in figure 9 
where the mean hydrogen density of 3 cm -3 is assumed in 
order to have a clump filling fraction larger than zero for 



3 As increasing the disc optical depth, the transmission rate grad- 
ually saturates at I„ in equation (21). If the layering parameter 
does not depend on the wavelength, the saturation value does 
not, either. Thus, we have a featureless and gray attenuation law. 



Figure 6. Dependence on the layering parameter for the diffuse 
stellar disc (§): (a) transmission rates and (b) attenuation laws 
normalized at visual band. In both panels, thick solid, dotted, and 
dashed curves are the composite transmission rates/attenuation 
laws with £ = 0.5, 1.0, and 2.0, respectively. In panel (a), thin 
solid, dotted, and dashed curves are transmission rates To. The 
transmission rate Ti is not affected by £. The Small Magellanic 
Cloud type dust is assumed. 
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Figure 7. Dependence on the inclination angle: (a) composite 
transmission rates and (b) composite attenuation laws normalized 
at visual band. The inclination angle is 0°, 29.4°, 42.1°, 61.1°, 
and 76.9° from top to bottom in both panels. The Milky Way 
type dust is assumed. 



Attenuation law through a clumpy ISM 9 




Figure 8. Dependence on the mean gas density: (a) composite 
transmission rates and (b) composite attenuation laws normalized 
at visual band. The hydrogen number density is 0.32, 0.82, 2.1, 
5.5, 14.4, and 37.3 cm -3 from top to bottom in both panels. The 
Small Magellanic Cloud type dust is assumed. 



Figure 10. Dependence on the dust-to-gas mass ratio (X>): (a) 
composite transmission rates and (b) composite attenuation laws 
normalized at visual band. In both panels, solid, dotted, and 
dashed curves are the cases with V = 10 -2 , 10" 2 ' 5 , and 10 -3 , 
respectively. The Small Magellanic Cloud type dust and the mean 
hydrogen density of 3 cm" 3 arc assumed. 
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Figure 9. Dependence on the equilibrium thermal pressure 
(p/te): (a) composite transmission rates and (b) composite at- 
tenuation laws normalized at visual band. In both panels, solid, 
dotted, and dashed curves are the cases with p/k B = 10 3 - 5 K 
cm" 3 , 10 4 K cm" 3 , and 10 3 K cm" 3 , respectively. The Milky 
Way type dust and the mean hydrogen density of 3 cm" 3 arc 
assumed. 



the p/k B = 10 3 K cm 3 case. If we adopted the mean den- 
sity of 1 cm" 3 , we would have nu = riH.wnm (see eq.l), and 
then, / c i = for the p/k B = 10 3 K cm" 3 case. The pres- 
sure affects on the attenuation law by two ways: the density 
contrast between clumps and the inter-clump medium, and 
the clump optical depth through the size (i.e. Jeans length). 
Approximately, the density contrast is proportional to p 0A 
found from equations (1) and (2), and the clump optical 
depth is proportional to p ' 5 found from equation (3). Thus, 
the effect of the pressure is not so large as shown in figure 9 
for the limited range of the value considered here. 

5.1.6 Dust-to-gas mass ratio 

Figure 10 shows the effect of the dust-to-gas mass ratio. Here 
we adopt the mean hydrogen density of 3 cm" 3 to avoid 
having a too small disc optical depth for the case of T> = 
1CP 3 , otherwise we have almost zero attenuation at visual 
band which is the normalization of the presentation. We find 
that a steeper attenuation law as T) decreases because the 
disc optical thickness decreases. Conversely, we can make a 
grayer attenuation law if V increases. However, it is unlikely 
that many galaxies have a much higher dust-to-gas ratio 
than that of the Milky Way. 

5.1.7 Global geometry 

We examine the effect of the global geometry. Suppose a 
sphere having a two-phase medium with an equilibrium ther- 
mal pressure and a mean gas density. Young stars are con- 
fined in clumps, and clumps and older stars distribute uni- 
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Figure 11. Comparison with a disc geometry and a spherical ge- 
ometry: (a) transmission rates and (b) attenuation laws normal- 
ized at visual band. The solid curves are the spherical case and the 
dotted curves are the disc case. The thin curves in the panel (a) 
are transmission rates To (upper) and T\ (lower). The thick curves 
in both panels are the composite transmission rates/attenuation 
laws. The Milky Way type dust is assumed. 

formly in the sphere. When we consider isotropic scatterings 
and assume that the scattered photons are also distributed 
uniformly in the sphere, the transmission rate for diffuse 
sources can be analytically expressed as equation (10) with 
the effective optical depth of K ffR and the effective albedo 
iv e ff. The radius of the sphere R is set to be equal to the 
dusty disc for the comparison. The clump optical depth 
and albedo are also given by equations (3) and (9). There- 
fore, we have the transmission rate of the clumpy sphere 
as 

T s P h = /yT s P h + (1 _ /y ) T Bph ; (25) 

where the transmission rate for the diffuse source is 

T sph = Posers h d , ui ef f ) , (26) 

and the transmission rate for the clumpy source is 

j*p _ p csc ( Tcl](iJcl )p csc ( Kcff /,, d)a;cff ) _ (27) 

Figure 11 shows the comparison of the disc and spher- 
ical geometries. We have assumed the standard set tabu- 
lated in table 1 in this comparison. For the disc geometry, 
we adopted £ = 0.5 since this choice is not so important 
as shown in section 5.1.2. We find very similar transmis- 
sion rates. Differences appear in longer wavelengths where 
the optical depth is small. This differences are caused by 
the scattering. In the face-on disc geometry, the path length 
of the photons scattered into the observer's (face-on) line 
of sight is shorter than the original length. On the other 
hand, the scattered path length is longer than the original 
one in the spherical geometry. Thus, the face-on disc is more 



Figure 12. Different dust property in clumps: (a) transmission 
rates and (b) attenuation laws normalized at visual band. The 
solid curves are the case that the extinction law of Ry = 5.5 is 
assumed only for dust in clumps, and the dotted curves are the 
case that the clump dust is the same as the diffuse dust (Ry = 
3.1). The thin curves in the panel (a) are transmission rates To 
(upper) and Ti (lower). The thick curves in both panels are the 
composite transmission rates/attenuation laws. The Milky Way 
type dust is assumed. 

transparent for the scattered photons than the spherical ge- 
ometry. This leads to the difference of the attenuation law 
normalized at visual band because of the different normal- 
ization. 

5.1.8 Grayer dust in clumps 

The observed extinction laws in the Milky Way and Mag- 
ellanic Clouds show a very large variety among lines of 
sight and are expected to depend on the environment 
(e.g., iFitzpatrickl Il998t iGordon et all 120031) . Especially, 
the dust property in star-forming regions (i.e. clumps) 
may be different from tha t in d iffuse medium, fndeed, 
ICardelli. Clayton, fc Mathisl lll989l) showed that the extinc- 
tion laws along some lines of sight toward dense clouds have 
a large value of Rv = Av/E(B — V); a grayer extinction 
law. For example, they found Rv = 5.5 toward the Orion 
nebula, in contrast, Rv = 3.1 for the average extinction 
law. Such a large Rv was attributed to the grain growth in 
dense medium. On the other hand, iRachford et alJ (120021) 
reported that there is no correlation between Rv and the 
molecular fraction along translucent lines of sight. Although 
a systematic difference of the dust property in dense medium 
(a grayer extinction law) is not completely established, the 
possible effect is worth examining. 

We assume that the dust in clumps is grayer than the 
diffuse dust; the extinction cross section (fed, a) in clumps is 
set to be the extinction law with Rv = 5.5 by Cardelli et 
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al. (1989). The extinction cross section for the inter-clump 
dust is the same as introduced in section 3.1. This cor- 
responds to the case of Rv ~ 3.1. The scattering albedo 
(uid.x) and the asymmetry parameter (gd,x) are not changed 
because we do not have enough informations about them 
for the case of Rv = 5.5. If the dust size of Rv = 5.5 is 
larger than that of Rv = 3.1, the adopted scattering albedo 
for Rv = 5.5 may be smaller than the real values in op- 
tical wavelengths. It leads a underestimation of the trans- 
mission rate Ti. However, the composite attenuation at the 
visual band is dominated by To which is determined by the 
Rv = 3.1 dust. Thus, the effect is not important. 

Figure 12 shows the case assuming the grayer dust of 
Rv = 5.5 in clumps as solid curves. For the dotted curves, 
the clump dust is the same as the diffuse dust which has 
always Rv = 3.1. From the top panel, we find that the 
transmission rate Ti of Rv = 5.5 becomes larger and grayer 
than that of Rv = 3.1. This makes the composite attenua- 
tion law grayer (bottom panel). However, it is still steeper 
than the extinction law (of Rv = 3.1) and the Calzetti law. 

5.2 Calzetti law 

The attenuation law obtained bv lCalzetti et, all l)l994h from 
observations of the UV bright starburst galaxies is grayer 
than the extinction laws of the Milky Way and the Magel- 
lanic Clouds. Here we try to reproduce the Calzetti law in 
the framework of this paper: the attenuation law through 
the two-phase ISM and the age-selective attenuation. 

Although we have obtained a steep attenuation law for 
normal disc galaxies in section 4, we can also make a grayer 
attenuation law as seen in section 5.1, for example, a larger 
mean density. Figure 13 shows a good agreement between 
the Calzetti law (dotted) with our composite attenuation 
law for the face-on case with a mean hydrogen density of 
5-10 cm -3 (solid). Active star-forming galaxies are likely to 
have a large amount of gas, and the gas den sity should ex- 
ceed a certain critical value llKennicuttll 989) . The hydrogen 
density of 5-10 cm -3 corresponds to 50-100 Mqpc~ 2 for the 
disc half height of 150 pc. Indeed, this gas column density 
greatly e xceeds the critica l density for the star formation 
found bv lKennicutl Jl98SD . 

There are two alternative ways to reproduce the 
Calzetti law by our model: a highly inclined disc and a 
system dominated by young stars. Even in the case of the 
mean hydrogen density of 1 cm -3 , the composite attenua- 
tion laws with an inclination angle of 75-80° (dashed curves 
in figure 13) also show a good agreement with the Calzetti 
law. However, this is unlikely because there is no systematic 
bias toward very hig h inclinations in the sample galaxies of 
ICalzetti et al.l <ll994l) . Since the galaxies really show an ac- 
tive star formation, the luminosity may be dominated by 
young stars. In this case, the attenuation law with f y = 1 
(i.e. attenuation law produced only by Ti) is suitable. The 
dot-dashed curves show such cases with a mean hydrogen 
density of 3-10 cm" 3 (face-on). Again, a larger density is 
preferred. 



5.3 Chariot and Fall prescription 

ICharlot fc Fail ( 2000) proposed a simple prescription of the 



attenuation law with the age-selective attenuation and ex- 
plained a lot of observational properties of the UV bright 
starburst galaxies. They adopted the age-selective attenua- 
tion by a similar way to us: the optical depth for old stars 
is only that by dust diffusely distributed in the ISM tism 
and the optical depth for young stars is the sum of tism and 
the optical depth by dust in the birth cloud tbc. Then, they 
simply assumed that the wavelength dependences of tism 
and tbc are the same: A -0 ' 7 . We examine this point. 

In figure 14, we show the normalized attenuation laws 
for Ti (dashed) and To (dotted) with several mean hydrogen 
densities, and compare them with the pow er-law attenua- 
tion w ith the index -0.7 (solid) proposed bv lCharlot fc Fall 
(2000). We find that (1) the slopes of the attenuation laws 
for the embedded young stars (Ti) and for the diffuse stars 
(To) are not the same. (2) The slopes decreases, especially 
for the diffuse stars, as the mean gas density increases. (3) 
The slopes for both stellar populations are similar to each 
other in a case with a mean hydrogen density of ~ 5 cm" 3 , 
in which the composite attenuation law is also similar to 
the Calzetti law. Therefore, the Chariot and Fall prescrip- 
tion and also the Calzetti law can be interpreted as a special 
case of a larger density (i.e. starburst) disc in the framework 
presented here. 



6 CONCLUSION 

The attenuation laws through the interstellar medium (ISM) 
with clumpy distributions of stars and dust was investigated 
in this paper. As an application of the multi-phase ISM 
model, we introduced a physical model of the dumpiness of 
the medium (i.e. dust distribution). In the model, the two 
parameters determining the dumpiness, the density contrast 
between clumps and the inter-clump medium and the vol- 
ume filling fraction of clumps, were translated into two phys- 
ical quantities in the ISM: the equilibrium thermal pressure 
and the mean gas density. The radiative transfer through the 
clumpy ISM was treated in a 1-D plain-parallel g eometry by 
the m ega-grain approximation developed bv lVarosi fc Dwekl 
(1999). For the dumpiness of the stellar distribution, we 
considered two stellar populations: older stars diffusely and 
smoothly distributed in the medium and younger stars em- 
bedded in clumps (i.e . birth clouds). T his provides the age- 
selective attenuation dSilva et alJll99g|) : young stars are sur- 
rounded by a denser dusty medium and suffer a stronger 
attenuation. The transmission rate of the entire radiation 
from the two stellar populations is a composite of the trans- 
mission rates of each individual population with a lumi- 
nosity weight. The stellar population dominating the radia- 
tion changes from older stars to younger stars as the wave- 
length decreases. With this fact, the age-selective attenu- 
ation makes the attenuation law steep; the composite at- 
tenuation rapidly increases from small attenuation for older 
stars at a long wavelength to large attenuation for younger 
stars at a short wavelength. Therefore, the attenuation law 
expected for normal disc galaxie s is much steeper than the 
Calzetti law (Calze tti et alJll994T) and the extinction laws in 
the Milky Way and the Magellanic Clouds. After an exten- 
sive survey of the parameter space, it was shown that the 
steep attenuation law is very robust for the parameter set 
expected in normal disc galaxies. Furthermore, the Calzetti 
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Figure 13. Comparison of the attenuation laws with the Calzetti law: (a) Milky Way type dust and (b) Small Magellanic Cloud type 
dust. The dotted curve is the Calzetti law. The solid curves are the composite attenuation laws for the face-on case with the mean 
hydrogen number densities of (a) 5.5 cm -3 and (b) 8.9 cm -3 . The dashed curves are the composite attenuation laws of the mean 
hydrogen density of 1 cm -3 for the inclination angle of (a) 76.9° and (b) 80.7°. The dot-dashed curves are the face-on attenuation law 
for Ti with the mean hydrogen densities of (a) 3.4 cm -3 and (b) 8.9 cm -3 . 




Figure 14. Comparison of attenuation laws with the prescription by Chariot & Fall (2000): (a) Milky Way type dust and (b) Small 
Magellanic Cloud type dust. The solid line is the power-law attenuation law by Chalrot & Fall (2000); oc A - 7 . The dotted and 
dashed curves are the attenuation laws corresponding transmission rates To and Ti , respectively. The mean hydrogen number densities 
are 0.51, 2.1, and 8.9 cm -3 from top to bottom for both types of curves and in both panels. 



law was reproduced in the framework presented here with 
a larger gas density than that in normal disc galaxies. Such 
a large gas density is very likely for the starburst galaxies. 
Finally, we examined the validity o f a simple pres c riptio n 
of the attenuation law proposed bv lCharlot fc Fall {2000). 
Then, we found that their prescription (and also the Calzetti 
law) is a special case with a larger density disc in our model. 

Interestingly, the UV color (or spectral slope) of nor- 
mal disc galaxies is systematically redder than that of UV 
bright starbursts at a fixed flux ratio o f the UV t o the in- 
frared (or fixed UV attenuation) (e.g., lBelJl2002l). A firs t 
result of GALEX also shows such trend l|Buat et aljEoOst) . 
One of the explanations of the trend is that normal disc 



galaxies have intrinsically r ed colors because of their qui- 
escent star-forming activity jKong et al.ll2004h . However, a 
steep attenuation law presented here can make galaxies red 
with small attenuation, and then, can be an alternative in- 
terpretation of the trend. Nowadays, GALEX is providing 
us with the UV data, including spectra, of a large number of 
normal disc galaxies. In the near future, GALEX will con- 
clude whether the attenuation law of normal disk galaxies is 
steep or not. 
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